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ABSTRACT 

The arched equilibrium path of a thin magnetic flux tube in a 
plane-stratified, nonmagnetic atmosphere is calculated for cases in which the 
flux tube contains a steady siphon flow. The large-scale mechanical 
equilibrium of the flux tube involves a balance among the magnetic buoyancy 
force, the net magnetic tension force due to the curvature of the flux-tube 
axis, and the inertial (centrifugal) force due to the siphon flow along curved 
streamlines. The ends of the flux tube are assumed to be pinned down by 
some other external force. Both isothermal and adiabatic siphon flows are 
considered for flux tubes in an isothermal external atmosphere. For the 
isothermal case, in the absence of a siphon flow the equilibrium path reduces 
to the static arch calculated by Parker (1975, 1979). The presence of a 
siphon flow causes the flux tube arch to bend more sharply, so that magnetic 
tension can overcome the additional straightening effect of the inertial force, 
and reduces the maximum width of the arch. The curvature of the arch 
increases as the siphon flow speed increases. For a critical siphon flow, with 
supercritical flow in the downstream leg, the arch is asymmetric, with greater 
curvature in the downstream leg of the arch. Adiabatic flows have 
qualitatively similar effects, except that adiabatic cooling reduces the buoyancy 
of the flux tube and thus leads to significantly wider arches. In some cases the 
cooling is strong enough to create negative buoyancy along sections of the 


flux tube, requiring upward curvature of the flux tube path along these 
sections and sometimes leading to unusual equilibrium paths of periodic, 
sinusoidal form. 


Subject headings: hydromagnetics — plasmas — Sun: atmospheric 
motions — Sun: magnetic fields — Sun: sunspots 
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I. INTRODUCTION 

In the two previous papers in this series (Thomas 1988, hereafter called 
Paper I; Montesinos and Thomas 1989, hereafter called Paper II) we 
explored the phenomenon of a siphon flow in a thin, isolated, arched 
magnetic flux tube in a stratified atmosphere under gravity. In those papers 
we ignored the problem of determining the equilibrium path of the flux tube in 
the atmosphere; instead, we assumed a shape for the flux-tube arch and then 
calculated the possible siphon flows along the flux tube. Now we relax that 
assumption and turn to the full problem of determining the exact equilibrium 
path of an arched flux tube containing a siphon flow. The siphon flow affects 
the equilibrium path of the flux tube by introducing inertial forces due to 
curvature of the flow streamlines and by altering the magnetic tension along 
the flux tube. 

The problem of determining the equilibrium path of a static, thin, isolated 
magnetic flux tube in a plane-stratified atmosphere was first solved by 
Parker (1 975, 1 979). He considered the cases where the external 
atmosphere is either isothermal or polytropic and, in both cases, he took the 
temperatures inside and outside the flux tube to be equal. He found 
equilibrium paths h(x) (height h as a function of horizontal distance x) in the 
form of a symmetric arch, provided the footpoints of the arch are anchored 
by some other force. In the equilibrium configuration, the magnetic buoyancy 
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of the flux tube is balanced everywhere along the tube by the net magnetic 
tension force due to the curvature of the flux tube axis. This curvature is 
everywhere concave downward (d 2 h/dx 2 < 0) in order to balance the 
upward buoyant force; thus, as we move away from the top of the arch, at 
some point the flux tube becomes vertical, beyond which the buoyant force 
cannot be balanced. This results in a limiting, maximum width of the flux tube 
arch, which works out to be equal to a few density scale heights of the external 
atmosphere. 

In the present paper we show that arched equilibrium paths also exist for 
thin flux tubes containing a siphon flow. As far as the equilibrium path of the 
flux tube is concerned, the most important new feature introduced by the flow 
is the inertial force (centrifugal force) due to flow along curved streamlines. 
The net magnetic tension force must balance this inertial force in addition to 
the buoyancy force, which requires that the arch have greater curvature 
than in the static case. This leads to a somewhat smaller maximum width of the 
arch. The siphon flow also affects the shape of the arch by altering the 
magnetic field strength along the tube through the Bernoulli effect (decreased 
internal gas pressure). A brief mention of these results has already been 
given elsewhere (Thomas 1989; Thomas and Montesinos 1989). 

Spruit (1981) has given a derivation of the equations describing the 
equilibrium path of a thin flux tube, using a method different from Parker's. 


6 


Browning and Priest (1984, 1 986) extended the work of Parker and Spruit by 
considering the effect of adding an ambient magnetic field to the external 
atmosphere. 

As in Papers I and II, the atmosphere outside the flux tube is assumed to be 
a plane-stratified, isothermal, perfect-gas atmosphere in hydrostatic 
equilibrium. The subscript e will be used to denote quantities associated with 
the external atmosphere. Letting T e0 denote the uniform exterior 
temperature, the distributions of pressure and density with height h outside 
the flux tube are given by 

P e = P e0 ex P(’ h/H ) - P e = P a o ex P(* h/H ) • O 1 - 1 ) 

where H = RT e0 /g is the constant scale height and p e0 and p e0 are the 

reference external pressure and density at h = 0. 

The main purpose of this paper is to show how the equilibrium path of a 
thin, isolated magnetic flux tube in a stratified atmosphere can be calculated in 
cases where the flux tube contains a siphon flow and to illustrate the effect of 
the siphon flow on the path of the flux tube. For this purpose, our assumption 
of a simple isothermal external atmosphere is sufficient. In a subsequent 
paper in this series we will adopt a more realistic models of the solar 
convection zone and photosphere as the external atmosphere in order to 
relate our siphon flow calculations more directly to intense magnetic flux 
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tubes on the Sun. 

Note: After this paper was submitted for publication, a paper by 
Degenhardt (1989) appeared in which he also solves for the equilibrium path 
of a thin, isolated magnetic flux tube containing a siphon flow, using essentially 
the same method we present here. Other than the basic method, however, 
there is little overlap in the two papers. Degenhardt concentrates on specific, 
dimensional cases in a more detailed plane-stratified model of the external 
solar atmosphere, whereas we have emphasized an exploration of the full 
range of parameter values for a simple external atmosphere (isothermal). 
Degenhardt considers only isothermal flows, whereas we consider both 
isothermal and adiabatic flows. Also, Degenhardt does not integrate the 
governing equations down to the point where the flux tube becomes vertical 
and thus does not determine the maximum width of the arch as we do here. 
Degenhardt shows that in the case of a temperature-stratified 
(nonisothermal) external atmosphere, a smooth transition from subcritical to 
supercritical flow can sometimes occur in the ascending part of the flux tube 
arch, whereas in an isothermal external atmosphere such a transition can 
occur only at the top of the arch. 

II. BASIC EQUATIONS 

Our derivation of the basic equations describing the equilibrium path of a 
thin flux tube containing a siphon flow generally follows the approach of 
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Spruit (1 981 ). The momentum equation for steady flow is 

R 2 1 

pv-Vv = ' V(P + ^ B- VB + pg . (2.1) 

The flux tube is in lateral pressure balance with the external atmosphere: 

p + £ - p. ■ < 2 - 2 > 

(The subscript e denotes quantities associated with the external atmosphere, 
while unsubscripted variables refer to the interior of the flux tube.) Then, 
assuming hydrostatic equilibrium of the external atmosphere, we can write 

V(P+ 8^ )== Vp e = p e g ’ 

and the momentum equation (2.1 ) can be written in the form 

pv V v = -^B'VB + (p-p e ) g . (2.3) 

We assume that the flux tube is thin, in the sense that its radius is much 
smaller than either its radius of curvature or the scale height H of the 
external atmosphere. Then we can consider the flow in the interior of the flux 
tube to be one-dimensional, i. e., locally parallel and uniform across the 
cross-section of the flux tube, so that all interior variables are functions only 
of the tangential coordinate s. Let o denote the unit vector along the axis of 
the flux tube (in the s-direction) and let v denote the unit vector normal to the 
axis of the flux tube, directed toward the center of curvature (see Fig. 1 ). 

Then we decompose the magnetic tension force into tangential and normal 
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components as follows: 

BVB = B £ (Bo) = (B§ )o + B 2 f = (l^f )a + B 2 k , (2.4) 

where k = da/ds = v/R is the radius of curvature vector, with R denoting the 
radius of curvature of the flux tube axis. Similarly, for the inertial force we 
can write 

2 

v-Vv = (T-Tr-)a + v 2 k . (2.5) 

v 2 ds 

The acceleration of gravity g can also be decomposed into tangential and 
normal components, in the form 

g = (-g sin6)o + (g cos6)Rk , (2.6) 

where 0 is the angle between the flux tube axis and the horizontal direction 
(see Fig. 1 ). Using expressions (2.4) - (2.6), we can write the a and k 
components of equation (2.3) as follows: 


p ^< y > - i ( ^)-^- p)9Sin9 = 0 ■ 
2 

R ( P y2 ' 4 ^" ) + (Pe - P)9 C °S9 = 0 • 


(2.7) 

( 2 . 8 ) 


The tangential component, equation (2.7), describes the siphon flow in the flux 
tube and, as we shall show below, is equivalent to the momentum equation 
used in Papers I and II. The normal component, equation (2.8), describes the 
large-scale mechanical equilibrium of the flux tube; it represents a balance 
among the inertial force pv 2 /R due to the curvature of the flow streamlines, 


the magnetic tension force - B 2 /4 tiR due to the curvature of the flux tube, and 
the normal component (p e - p)g cos0 of the buoyancy force. 

If we let h(x) denote the height of the flux tube axis (see Fig. 1 ), then we can 
write 


d_h 

dx 


R 


= tane 


d 2 h 

dx 2 


i de 


2 fl dx ’ 


d 2 h 

dx 2 


1+ 0 2 


-i3/2 


cos 0 


-COS0^- . 
dx 


With these relations, we can formulate the problem of finding the equilibrium 
path h(x) of the flux tube, based on equation (2.8), as a pair of coupled 


first-order differential equations for 0(x) and h(x): 


zr- = tan0 , (2.9) 

dx 

(^--pv 2 )^-= -(p e -p)g . (2.10) 

These two equations together with the pair of first-order differential 
equations for the velocity v and cross-sectional area A of the siphon flow 
(derived in Papers I and II) and the condition of magnetic flux conservation, AB 
= constant, provide a set of equations for determining simultaneously the 
siphon flow inside an arched flux tube and the equilibrium path h(x) of the flux 
tube. 


The tangential component of the steady momentum equation can be written 


in the form 


1 1 


,dv 


dp 


.dh 


^ v ds ds ^ds 


( 2 . 11 ) 


noting that sin0 = dh/ds. Equation (2.1 1), which is the form of the momentum 
equation used to describe the one-dimensional siphon flows in Papers I and 
II, follows directly from equation (2.1), and also follows from equation (2.7) if 
we again make use of the condition of lateral pressure balance, equation 
(2.2), and the hydrostatic equation for the external atmosphere, Vp e = p e g. To 
complete the description of the one-dimensional steady siphon flow within the 


flux tube, we also need the equation of mass conservation, 

(pvA) = 0 , pvA = Q = constant, 

the equation of magnetic flux conservation, 

-7- (AB) = 0 , AB = <J> = constant , 

ds ' ' 


the equation of state, 


p = pRT , 


and the appropriate energy equation, 


isothermal flow: T = T c 




adiabatic flow: 




( 2 . 12 ) 


(2.13) 


(2.14) 


(2.15) 


Equations (2.1 1 ) - (2.1 5), together with equation (2.2), constitute the complete 
set of equations describing the siphon flow, as used in Papers I and II. As 
before, we are using the isothermal and adiabatic cases to represent 
extremes between which the behavior of a flow with a more realistic energy 


equation (including radiative exchange) will lie. 


It is convenient to write the basic equations in terms of dimensionless 
quantities. We adopt the same scaling as in Papers I and II, except that here 
we use the scale height H of the external atmosphere as the length scale for 
both the height h and the horizontal distance x. In Papers I and II, the 
horizontal distance x was scaled by the distance L between the footpoints of 
the arch, and L was left undetermined. Here it is more convenient to use the 
same scaling for h and x, so that dimensionless plots of h(x) show the actual 
shape of the flux tube arch. The actual value of the distance L between 
footpoints is determined by the equilibrium solution for h(x). Also, here we 
take the origin x = 0 to be at the point where the arch reaches its maximum 
height (see Fig. 1), whereas in Papers I and II the origin was taken to be at the 
left footpoint of the arch (where h = 0). 

Using H as the length scale and the exterior isothermal sound speed c, = 

(RT^) 1 / 2 as the velocity scale, we define the following dimensionless quantities 
(denoted by an overbar): 



where p 0 , p 0 , T 0 , A 0 , B 0 , and v 0 are the reference values (at h = 0 at the 

upstream footpoint) of the pressure p, density p, temperature T, 
cross-sectional area A, magnetic field strength B, and velocity v inside the flux 


tube. The AIfv6n speed a, the "tube speed" c t , and the characteristic speed c 1 
are defined by 


a 


2 




^2 
c a 

2 2 
c + a 



(2.17) 


as in Papers I and II. The sound speed c inside the flux tube is taken to be the 
isothermal sound speed (RT) 1 / 2 for isothermal flow and the adiabatic sound 
speed (yRT) 1 / 2 for adiabatic flow. We also define the "plasma beta" based on 
the external pressure as follows: 


m - 


8K Pe (h) 

B 2 (h) 


8rcp 


p 0 - p(°) = , 


eO 


(2.18) 


where p o (called simply (3 in Papers I and II) is the particular value of (3 at h = 
0. In terms of p, the lateral pressure balance condition (2.2) may be 


expressed as 


P P-1 

Pe " P 


(2.19) 


III. ISOTHERMAL FLOWS 


In this section we consider isothermal flows, for which the temperature 
inside the flux tube is everywhere equal to the uniform external temperature 
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i. e., T = T e0 = constant. In this case the condition of lateral pressure balance, 

equation (2.2), allows us to express the buoyancy term on the right hand side 
of equation (2.10) as 

•(p e -p)g = * RT (p e ■ p) * 

eO 

Then, in terms of the nondimensional variables defined in equations (2.16) 
and (2.17), the pair of equations (2.9) and (2.10) can be written in the form 

= tan0 , (3.2) 

dx 

( 3 . 3 ) 


As shown in Paper I, the set of equations (2.2) and (2.11) - (2.15) in the case 
of isothermal flow lead to the velocity-height relation 


dv _ -i 2-2 - dh 

dx - v \2-9 flS 0 A exp(-h) + 1]J d x 


(3.4) 


and the area-height relation 


Equations (3.2) - (3.5), together with the simple condition of magnetic flux 
concentration, AB = 1 , constitute a coupled set of nonlinear equations for 
determining simultaneously the siphon flow variables v(x) and A(x) and the 
equilibrium path h(x) of the flux tube. However, computationally the problem 
is actually decoupled in the sense that equations (3.4) and (3.5) can be written 
with height h as the independent variable: 


{ _J 

[2-v 2 [P 0 A 2 exp(-h) + 1]J ’ 

- , 1 - v 2 p.A 2 exp(- h) 

2 :° 2 ~ 
dh ^ 2 — v [P 0 A exp(- h) + 

Thus, we can solve for the flow variables v(h) and A(h) using equations (3.6) 
and (3.7), subject to imposed values at some height h. The other flow 
variables p, p, T, and B are then given as functions of h by using equations 
(2.12) - (2.15). We can then solve equations (3.2) and (3.3) for the equilibrium 
path h(x) of the flux tube using the known quantities v(h), B(h), and p(h). (We 
should point out that this decoupling of the problem is a direct consequence of 



the thin flux tube approximation; in a thick flux tube, the flow will not be 
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one-dimensional and the problem of determining the equilibrium path and the 
flow will be fully coupled.) 

In order to solve for both the siphon flow and the equilibrium path, we 
proceed as follows. We choose a value of p 0 and specify the nondimensional 

height a of the arch above the reference level h = 0. Then we can determine 
the siphon flow variables by integrating equations (3.6) and (3.7) beginning at 

h = 0, where we have the specified initial values v(0) = v 0 , A(0) = 1 . We use a 

fourth-order Runge-Kutta method with variable mesh spacing (Press etal. 
1986) for the numerical integration. Once the siphon flow is determined up to 
the top of the arch and we know v(h) and B(h), equations (3.2) and (3.3) can be 
integrated starting from the top of the arch, where we also have the known 
initial conditions h = a , 6 = 0. Also, we can take the values of v and A at the top 
of the arch determined by the numerical integration of equations (3.6) and 
(3.7) and use them as initial values for the equations ((3.4) and (3.5). In this 
way, the full set of four equations (3.2) - (3.5), with the condition AB = 1 , can 
be integrated simultaneously, beginning at the top of the arch (at x = 0), in 
order to determine h(x), v(x), A(x), and B(x), and hence all other flow 
variables. The integration can be carried beyond the point where h = 0, 
since the choice of the reference level h = 0 is arbitrary. The integration can 
proceed up to the point where the flux tube becomes vertical, i. e., where 


dh/dx becomes infinite. This yields a full solution of the problem for the 
parameter values a, p o , and v 0 . We are free to choose any other height 

(below the top of the arch) as the reference level h = 0 for this solution, in 
which case the solution corresponds to a different arch height a and different 

values of p 0 = p(0) and v 0 = v(0). In fact, we could just begin by specifying the 

values of v, A, and p at the top of the arch and integrate equations (3.2) - (3.5) 
from there, but the method we have chosen is more consistent with our 
approach in Papers I and II. 

In the absence of a siphon flow (v = 0), the pair of equations (3.2) and (3.3) 
for determining the path h(x) of the flux tube reduce to 


4=- = tane , £--l 

dx dx 2 


(3.8) 


These equations are equivalent to those solved by Parker (1975, 1979) and 


have the analytical solution 


tan 2 (x/2) = exp(a - h) - 1 (3.9) 

for the case where h = a when x = 0. The resulting arch shape h(x) is shown 
in Figure 2a ( dot-dash curve), where we have taken a = 1 .0. Note that h 

as x ->• ± ji, so that the arch has a limiting maximum width L max = 2 ttH. Note 
also that in the static isothermal case, the equilibrium path is independent of 
the value of P 0 , i. e., independent of the magnetic field strength of the tube. 


This is due to the fact that each of the two competing forces, the magnetic 
buoyancy and the net magnetic tension, scale as B 2 . With a siphon flow, 
however, the added inertial force does not scale simply as B 2 and the 
equilibrium path does depend on |3 0 (see equations (3.2) and (3.3)). 

In addition to the static case, Figure 2a shows examples of the computed 
equilibrium path of the flux tube for two subcritical siphon flows and the critical 
siphon flow (with both subcritical and supercritical downstream branches) 
for the parameter values a = 1 .0 and p 0 = 3.0. The corresponding 

distributions of velocity v(x), cross-sectional area A(x), and plasma beta P(x) 
for these flows are also shown in Figure 2. Note in Figure 2c that the critical 
flow produces bulge points (points of local maximum cross-sectional area) in 
the flux tube (see Paper I for a discussion of these bulge points). In Figure 2d 
we see that the flow reduces the value of P(x) (i. e., increases the magnetic 
field strength), compared to the static case, everywhere above h = 0. This is a 
consequence of the decreased internal pressure associated with the siphon 
flow (the Bernoulli effect). Here we have taken the value of p to be 3.0 at h = 0 
in all four cases (the static case and the three cases with siphon flows). To see 
the effect of the siphon flow on an individual flux tube of fixed total magnetic flux, 
it is more appropriate to take the limiting value of p as h -» - « (where v -» 0) 
to be the same in all cases; then the siphon flow will reduce the value of p 


below that of the static tube everywhere along the tube. 

Comparing the different curves of h(x) in Figure 2a, we see that for fixed 

P 0 increasing the siphon flow speed causes the flux tube to bend more 

sharply, in order that the net magnetic tension force can balance the sum of 
the inertial and buoyancy forces. The siphon flow thus reduces the maximum 
width of the arch below the value for the static arch. For example, in Figure 
2 a the critical flow with its subcritical downstream branch reduces the 
maximum width of the arch to L max = 5.31 H, compared to L max = 2 tiH = 6.28H 

in the static case. For the case of the critical flow with its supercritical 
downstream branch, the equilibrium path is asymmetric, with greater 
curvature in the supercritical downstream leg of the arch and a further 
reduction of the maximum width. Of course, the flow will not actually follow the 
supercritical branch all the way to the downstream footpoint; at some point in 
the downstream branch the flow will decelerate to subcritical speed through 
a standing tube shock and the curvature will be reduced beyond that point. 

We intend to present computations of critical flows with standing tube shocks in 
a subsequent paper in this series. 

The purely subcritical steady flows we present here (and in the next 
section) are all symmetric about the top of the arch and, in particular, the gas 
pressure is the same at both footpoints. If we were to add the effects of 
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viscosity, these subcritical flows would no longer be symmetric; the steady 
flow would require a smaller pressure at the downstream footpoint in order 
to balance the viscous stresses. There is no ambiguity about the symmetric, 
steady subcritical flows in the inviscid case, however, as suggested by 
Degenhardt (1989). These flows simply represent the steady-state limit of a 
transient flow that arises as a consequence of an initial pressure difference 
between the two footpoints, with the flow directed toward the lower-pressure 
footpoint. The critical inviscid flows with supercritical flow in the downstream 
leg of the arch are inherently asymmetric, although they too could have equal 
pressures at the two footpoints with a tube shock of suitable strength in the 
downstream leg. 

In Figure 3 we show the equilibrium path h(x) , velocity v(x), and 
cross-sectional area A(x) for critical flows with three different values of (3 0 

(2.0, 3.0, and 6.0), all with a = 1 .0, along with the static case for comparison. 
For each critical flow both the subcritical and supercritical downstream 

branches are shown. Note that for increasing (3 0 , i. e., for decreasing 

magnetic field strength, the critical siphon flow has greater effect on the 
equilibrium path of the arch (produces greater curvature), even though the 

flow speeds are lower for increasing p o . In Figure 3a, the maximum width of 
the arch, based on the subcritical downstream branch of the flow, is L max = 
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5.68H, 5.31 H, and 4.67H for the critical flows with p o = 2.0, 3.0, and 6.0 
respectively, which can be compared to the value L max = 6.28H for the static 

case. The width of the arch is even further reduced for supercritical 
downstream branches of the flows, but again in these cases the flow must 
decelerate abruptly to subcritical speed somewhere in the downstream 
branch across a standing tube shock, which we have not computed here. 


IV. ADIABATIC FLOWS 

Our treatment of adiabatic siphon flows follows the notation of Paper II. We 
use the second form of the energy equation (2.15), and we allow the 

temperatures T 0 and T e0 inside and outside the flux tube at the upstream 

footpoint to be different, letting x = (To/T e0 ) denote their ratio. Equations (2.9) 
and (2.1 0) for the equilibrium path can be written in the form 


dh 

dx 


“ = tanQ 


r Or —2 

KttVJB -pv |de 

L P ° J d * 


P 0 T 

fe\)exp(-h) ■ p 
p o 


(4.1) 

(4.2) 


where we have also made use of equation (1.1). As shown in Paper II, the set 


of equations (2.2) and (2.1 1 ) - (2.1 5) lead to the coupled equations 
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dv 

dh 


r-c 2 

L r (c t 2 -v\ 


dA 

dh 



(4.3) 


(4.4) 


for the velocity and cross-sectional area, with 


r= 


f 2y(c 2 + a 2 )] 

1-2 _2 I 

I 2c +ya | 


-2 _ Y -1 

c =typ 


-2 

a 


2x 

(P 0 -1)pA 2 


-2 c 2 a 2 

c t -im 
c +a 


-2 —2 -2 
_ 2 (2c + ya -2y)c 

C 1 = ~2 Ti 

2c + ya 




(4.5) 


(Note: Equation (3.8) of Paper II, which corresponds to equation (4.3) above, 
contains a misprint, and equation (4.3) above is the correct form of the 
velocity equation. The correct form of the equation was used for all of the 
calculations in Paper II and the results were not affected by the misprint.) 
Equations (4.3) and (4.4), together with the auxilliary relations (4.5), constitute 
a closed set of equations for the siphon flow as a function of height h; they 
replace equations (3.6) and (3.7) used in the case of isothermal flow and can 
be integrated to determine the flow conditions at the top of the arch. If we 


rewrite them with x as the independent variable, equations (4.3) and (4.4) 
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together with equations (4.1 ) and (4.2) form a closed set of equations for 
determining simultaneously the siphon flow and the equilibrium path h(x), by 
integrating from the known initial conditions at the top of the arch. 

For the adiabatic flows considered here, the limiting static case (for v = 0) 
corresponds to a flux tube with static adiabatic internal stratification, 
embedded in the isothermal external atmosphere. The internal pressure, 
density, and temperature in the static case are given by 

y i 

r i 7-1 r i 7_1 r i 

p.ll-(^KI , , T = h-(^)hl , (4.6) 

and the cross-sectional area and magnetic field strength are given by 

.1 1 

r - _f 2 _ r - _i 2 

A = |P 0 ex P('h) - (P 0 - 1 ) PJ , B='P 0 ex P(-h) -(P 0 -1)Pj (4.7) 

(see Paper II). Unlike the static isothermal tube, the static adiabatic tube has a 

limiting maximum height h| im = yr/(Y- 1) where the internal pressure, density, 

and temperature go to zero; we are thus limited to considering arch heights 

a < h |im (= 2.5 for t = 1 .0 and y = 5/3). 

In a sense, the static adiabatic case is unrealistic for solar applications. It is 
reasonable to consider a siphon flow to be nearly adiabatic if the flow speed is 
so great that there isn't enough time for the internal gas to be heated to 
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thermal equilibrium by thermal conduction (or radiation) from the 
surroundings. If the radiative exchange time for the tube is much longer than 
the time it takes the flow to traverse the arch, then the flow will be nearly 
adiabatic and there will be significant cooling of the internal gas as it expands 
in the rising part of the arch. However, in the case of slow flows or, in the 
extreme, a static tube, heat exchange with the surroundings will be important 
and the primary justification for treating the flow as adiabatic does not apply. 
Nevertheless, the static adiabatic tube defined by equations (4.6) and (4.7) is a 
physically possible state if the flux tube were truly insulated from its 
surroundings and is also the correct mathematical limit of the adiabatic flow 
solutions for v 0; as such, it needs to be understood. 

In the adiabatic case (either static tube or with a siphon flow) the 
temperature inside the flux tube can be greater than or less than the 
temperature of the surrounding atmosphere, and this temperature 
difference can enhance or reduce the buoyancy of the flux tube. For 

example, in a static adiabatic tube with x = 1 .0 the temperature inside the tube 
equals the external temperature only at the level h = 0; as we go above h = 0 
the tube is progressively cooler than the surroundings, and as we go below h = 

0 the tube is progressively hotter than the surroundings. In some cases the tube 
may reach a height (call it h') where the density increase due to adiabatic 
cooling completely offsets the density decrease due to the reduced internal 
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pressure (the magnetic buoyancy effect) and the tube becomes heavier than the 
surroundings above that height; then, the curvature of the tube must switch to 
being concave upward above the height h\ and the point h = h' is an inflection 
point of the curve h(x). The relative effect of the adiabatic cooling is greater for 

larger values of p 0 because the magnetic buoyancy decreases with increasing 

P°- 

Figure 4 shows examples of the equilibrium path of a static adiabatic flux 
tube in an isothermal external atmosphere, computed from equations (4.1) and 

(4.2) with the appropriate static distributions of p and B given in equations 

(4.6) and (4.7). In Figure 4a we have plotted the paths h(x) of tubes with 

several values of p 0 , but all with a = 1.0 (same maximum height) and x = 1.0 

(internal temperature equal to external temperature at h = 0). Each flux tube is 
cooler than the surroundings everywhere above h = 0 and hotter than the 

surroundings everywhere below h = 0. As p 0 increases, the relative effect of the 

adiabatic internal stratification becomes more pronounced and the density in 
the upper part of the arch becomes more nearly equal to the external density; 
this is manifested in a reduced curvature of the arch, and the top part of the arch 

becomes progressively flatter for increasing p 0 . Eventually p 0 reaches a value 
(somewhere between 4.5 and 5.0) where the density at the top of the arch is 
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exactly equal to the external density and the tube is neutrally buoyant there. 

Then, with a further increase in p 0 , the tube becomes heavier than the 

surroundings at the top of the arch and the arch must become concave upward 
there. The equilibrium solution then abruptly switches over to a curious, 

sinusoidal-like path, such as the paths shown for |3 0 = 5.0 and 6.0 in Figure 4. 

In these cases, the negative buoyancy of the lower parts of the sinusoidal path 
(concave upward) is balanced by the positive buoyancy of the upper parts of the 
sinusoidal path (concave downward). There is positive buoyancy in the upper 
parts of the tube in these cases because, whereas the internal density is 
decreasing algebraically with height [cf. eqs (4.6)], the external density is 
decreasing exponentially [cf. eqs. (1 .1 )], so eventually a height is reached 
where the internal density is again lower than the external density and the 

buoyancy is again positive. Note that the sinusoidal flux tubes with p 0 = 5.0 and 
6.0 in Figure 4 are everywhere cooler than their surroundings, since they don’t 

extend down to h = 0 where we have assumed t = 1 .0. Also note that for these 

two flux tubes the values of p 0 » P(0) are somewhat irrelevant because the 

tubes do not extend down to h = 0; unlike the static isothermal case, p is not 
constant with height in a static adiabatic tube, but instead varies as 


which follows from equations (1.1), (4.6), and (4.7). 

In Figure 4 b we show the equilibrium paths of static adiabatic tubes with 

greater maximum height (a = 2.0) than the tubes in Figure 4a. The behavior is 

similar, except that here there is an intermediate case (for (3 0 = 6.0) with 

sections of upward concavity (negative buoyancy) on either side of the arch 
but with downward concavity (positive buoyancy) at the top of the arch. This 
situation is possible because the curves of internal density (algebraic) and 
external density (exponential) as functions of height will cross at two different 
heights for some range of values of the reference densities at h = 0, provided 
the arch is high enough. 

Figure 5 shows computed examples of the equilibrium path of a flux tube 
containing an adiabatic siphon flow, together with the velocity and cross 

sectional area along the tube, for the parameter values a = 1 .0, |3 0 = 3.0, and x 
= 1 .0. Shown here are the static case (for v 0 = 0), two examples of subcritical 

flow (for v 0 = 0.1 5 and 0.25), and the critical flow (for v 0 = 0.31 6207) with its 

subcritical and supercritical downstream branches. As in the isothermal 
case, we see that increasing the siphon flow speed increases the curvature of 


the flux tube and reduces the overall width of the arch. Here the overall width 
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of the arch is reduced from L = 8.94H in the static case to L = 7.1 1 H in the 
case of the symmetric critical flow with its subcritical downstream branch (the 
long-dashed curve). Comparing Figure 5a with Figure 2a, we see that the 
adiabatic case produces wider arches than the isothermal case. The 
adiabatic cooling of the gas leads to greater internal densities (compared to 
the isothermal case) and hence a weaker buoyancy force, which in turn 
requires less curvature of the flux tube axis in order to establish mechanical 
equilibrium. 

As expected, the critical adiabatic flow exhibits bulge points (Fig. 5c). The 
somewhat curious fact that the siphon flow speed (Fig. 5b) reaches a minimum 
and then begins to increase rapidly as we go down either leg of the arch is a 
consequence of the fact that the tube gets progressively hotter as we go lower 
down (because of the adiabatic assumption). Thus, the "local" value of the 
temperature ratio x (as opposed to the fixed value x = 1 .0 at h = 0) increases 

with depth and eventually exceeds the critical value x 2 s [y(p-l )+2]/yP and 

there the flow behaves as in the "superheated " tubes discussed in Paper II 
(see the discussion in section IV and Fig. 9a of Paper II). Of course, it is 
unrealistic to apply the adiabatic assumption all the way down the arch as h 
-oo, so the large increase in velocity there is artificial. 

In Figure 6 we show flux tubes with adiabatic siphon flows similar to those in 
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Figure 5, except with (3 0 = 6.0 instead of 3.0. For the critical flow ( with v 0 = 

0.228438) and for the purely subcritical flow with v 0 = 0.2, there is an 
equilibrium path in the typical form of an arch of limited width. However, if the 
initial velocity is further reduced to v 0 = 0.15, there is no longer an arched 

equilibrium path extending downward without limit. Instead, the integration of 
the basic equations away from the center x = 0 gives a sinusoidal-shaped 

equilibrium path, as does the static case ( v 0 = 0). Here we have a case of a 
flux tube that requires a siphon flow of sufficient speed to allow a simple 
arched equilibrium path; if the flow is too slow, then there is no simple arched 
equilibrium path, but there is an alternative sinusoidal equilibrium path with 
the same flow speed at the central point x = 0. For the sinusoidal equilibrium 

paths, which do not extend down to h = 0, the value of v 0 is somewhat 

irrelevant; it is perhaps better to think of the examples in Figure 6 as a 
sequence of flows with different specified values of the flow speed v at the 
central point x = 0. 


V. CONCLUSIONS 

We have shown how to calculate the equilibrium path of a thin magnetic flux 
tube in a stratified, nonmagnetic atmosphere when the flux tube contains a 
steady siphon flow. The equilibrium path is affected by the siphon flow, but the 
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problem is decoupled in the thin flux tube approximation because the siphon 
flow variables can be determined as functions of height h irrespective of the 
actual equilibrium path h(x). The large-scale mechanical equilibrium of the 
thin flux tube involves a balance among the magnetic buoyancy force, the net 
magnetic tension force, and the inertial force (centrifugal force) due to the 
siphon flow along a curved path. 

In general, the equilibrium path of a static thin flux tube in an infinite 
stratified atmosphere takes the form of a symmetric arch of finite width, with 
the flux tube becoming vertical at either end of the arch. A siphon flow within 
the flux tube increases the curvature of the arched equilibrium path, in order 
that the net magnetic tension force can balance the inertial force of the flow, 
which tries to straighten the flux tube. Thus, a siphon flow reduces the width of 
the arched equilibrium path, with faster flows producing narrower arches. 

The effect of the siphon flow on the equilibrium path is generally greater for 
flux tubes of weaker magnetic field strength. 

We have shown examples of the equilibrium for both isothermal and 
adiabatic siphon flows in thin flux tubes in an isothermal external atmosphere. 
For isothermal flows, the magnetic buoyancy is always positive, and the 
equilibrium path h(x) always has downward curvature (d^h/dx 2 < 0) and the 
equilibrium path must be an arch of finite width. For adiabatic flows, adiabatic 
cooling reduces the buoyancy of the flux tube and produces arches of 
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greater width than those with isothermal flows. In some cases adiabatic 
cooling more than offsets the magnetic buoyancy effect and causes sections of 
the flux tube to be heavier than its surroundings (negatively buoyant) and thus 
to have upward curvature (d 2 h/dx 2 > 0); this leads to the possibility of 
sinusoidal equilibrium paths of infinite horizontal extent. 

The isothermal and adiabatic siphon flows computed in Papers I and II 
were based on an assumed equilibrium path in the form of a parabolic arch 
of unspecified width. It turns out that the parabolic paths of the examples in 
Papers I and II are very close to the true equilibrium path computed by the 
methods of the present paper when we match the width of the arch in the two 
cases. Thus, the examples in Papers I and II are fairly accurate as they stand 
(within a few percent), with their assumed parabolic paths. However, since 
the siphon flow variables actually depend only on the height h and not on its 
distribution h(x), the flows computed in Papers I and II are precisely accurate 
if, instead of associating them with the assumed parabolic paths, we associate 
them with the exact equilibrium paths computed by the methods of the present 
paper. This can be thought of as mapping the siphon flow along the parabolic 
arch onto the true equilibrium arch by associating points of the same height h. 

The motivation for our study of siphon flows in isolated magnetic flux tubes 
is the possible applications to magnetic structures in the solar atmosphere, 
including intense photospheric flux tubes and sunspot penumbral filaments. 
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We have discussed these applications in some detail in Papers I and II. Here 
we note that in general the arched equilibrium paths of an isolated magnetic 
flux tube are fairly narrow, with footpoints separated horizontally by only 
some five to ten density scale heights (on the order of a thousand kilometers 
on the surface of the Sun). This short span of the flux-tube arch poses a 
problem for applications to penumbral filaments and intense photospheric 
flux tubes. Penumbral filaments are observed to extend horizontally over 
several thousand kilometers, for example. We have seen that adiabatic flows 
produce generally wider arches than isothermal flows, so to the extent that 
more realistic radiative flows allow some adiabatic cooling, somewhat wider 
arches may occur. 

However, a more likely mechanism for allowing the wider arches 
observed on the Sun has to do with other "canopy" magnetic fields at higher 
levels in the external atmosphere. Although we have considered the 
equilibrium path of a single, isolated flux tube without regard to any other 
magnetic fields in the external atmosphere, on the Sun the atmosphere is 
generally filled with magnetic field above some height in the low 
chromosphere. In the case of a penumbral filament this field is the general 
spreading canopy of the overall sunspot magnetic field, and in the case of an 
intense flux tube in the quiet photosphere this field is the general canopy 
magnetic field covering the inter-network regions of the quiet Sun (Giovanelli 
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1 980). We can imagine an isolated magnetic flux tube arching up to a point 
where it comes into contact with the overlying canopy field and then running 
horizontally for some distance along the base of the canopy before arching 
down again to the opposite footpoint (see also Thomas and Montesinos 1989). 
Along the section of the tube in contact with the canopy field, the mechanical 
equilibrium is established by a balance between the upward buoyancy of the 
flux tube and the downward magnetic pressure force due to the overlying 
canopy field. (Put another way, the isolated flux tube becomes part of the 
canopy field over this portion of the tube.) In this way, equilibrium arches 
(with flat tops) of much greater horizontal extent can occur on the Sun. The 
flow variables (velocity, cross-sectional area, etc.) will be quite uniform along 
the horizontal top of the arch in contact with the canopy, but in the isolated 
ascending and descending parts of the arch they will vary in the manner we 
have computed in this series of papers. Simple theoretical models of such 
extended equilibrium paths can be envisioned. 
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Figure Captions 

Fig. 1 . Basic geometry of an arched thin magnetic flux tube, showing the 
coordinate system and the unit tangential and normal vectors a and v. 

Fig. 2. Plots of (a) the equilibrium path h(x), ( b ) the velocity v(x), (c) the 

cross-sectional area A(x), and (d) the plasma beta P(x) for a thin magnetic 
flux tube containing isothermal siphon flows, with p 0 = 3.0 and a = 1 .0. 
Shown here are the static isothermal case ( dot-dash curves), which is 
independent of p 0 , and three examples of siphon flows: two subcritical 
flows ( solid curves) and the critical flow ( dashed curves) with both the 
subcritical ( long dash) and supercritical {short dash) downstream 
branches. The specified values of v 0 (the velocity at h = 0) for the three 

siphon flows are 0.1 5, 0.25, and v 0c = 0.303440. The corresponding 
maximum widths of the arch are L max = 6.05H and 5.69H for the purely 
subcritical flows and L max = 5.31 H for the critical flow (with its subcritical 
downstream branch), which can be compared with the value L max = 2 jcH = 
6.28H for the static case. 

Fig. 3. Plots of (a) the equilibrium path h(x) and (b) the velocity v(x) and 
cross-sectional area A(x) of the flux tube for critical isothermal siphon 
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flows with a= 1.0 and three different values of p 0 : p 0 = 2.0 {solid curves), 

P 0 = 3.0 ( dashed curves), and (3 0 = 6.0 (dotted curves). In each case both 

the subcritical and supercritical downstream branches are shown. Also 
shown for comparison is the static case (dot-dashed curve), which is 

independent of p 0 . The critical initial velocities for these flows are v 0c = 
0.339169, 0.303440, and 0.244038 and the maximum widths of the 
arches (for the subcritical downstream branches) are L max = 5.68H, 

5.31 H, and 4.67H for p o = 2.0, 3.0, and 6.0, respectively. 

Fig. 4. (a) The equilibrium paths of static adiabatic flux tubes in an isothermal 
external atmosphere, with a = 1 .0, x = 1 .0 (i.e., T 0 = T e0 ), and several 

different values of p 0 . For p 0 somewhere between 4.5 and 5.0, the flux 
tube becomes negatively buoyant at the top of the arch, and the equilibrium 
path for larger p 0 switches over to the sinusoidal shape represented 

here by the cases p 0 = 5.0 and 6.0. The maximum width of the arch is L max 
= 7.32H, 8.94H, 12.27H, and 17.1 8H for P 0 = 2.0, 3.0, 4.0, and 4.5, 
respectively, (b) Similar to Fig. 4a, but with a = 2.0. Here the case p 0 = 6.0 
has sections with upward concavity due to negative buoyancy on either 
side of the arch, and in the case p 0 = 7.0 the equilibrium has switched over 
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to a sinusoidal path. The maximum width of the arch is L max = 8.1 7H, 9.70H, 
and 15.69H for p 0 = 2.0, 3.0, and 6.0, respectively. 

Fig. 5. Plots of (a) the equilibrium path h(x), (b) the velocity v(x), (c) the 
cross-sectional area A(x), and (d) the plasma beta |3(x) for cases of 
adiabatic siphon flows with a = 1 .0, p 0 = 3.0, and x = 1 .0. Shown here are 

the static adiabatic case {dot-dashed curves), two examples of subcritical 
flows {solid curves), and the critical flow {dashed curves) with both its 
subcritical and supercritical downstream branches. The specified values 
of v 0 (the velocity at h = 0) for the three siphon flows are 0.15, 0.25, and v^ 

= 0.316207. The corresponding maximum widths of the arch are L max = 

8.60H and 7.94H for the two purely subcritical flows and L max = 7.1 1 H for 
the critical flow (with the subcritical downstream branch), which can be 
compared with the value L max = 8.94H for the static case. 

Fig. 6. Similar to Figures 5a, b, except with p 0 = 6.0. The critical flow ( dashed 

curves, with v 0c = 0.228438) and the faster subcritical flow (with v 0 = 0.20) 
allow a simple arched equilibrium path of finite width. The maximum width 
of the arch is L max = 1 8.30H for v 0 = 0.20 and L max = 1 1 .1 8H for the critical 
flow (with the subcritical downstream branch). The slower subcritical flow 
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(with v 0 = 0.15) and the static case {dot-dashed curves) do not allow a 

simple arched equilibrium path; instead, there is a sinusoidal-shaped 
equilibrium path with the appropriate flow speed at the center x = 0. The 
four cases here may be thought of as a sequence of flows with increasing 
flow speeds at the central point x = 0 (v = 0.0, 0.319, 0.481 , and 0.730, 
respectively, at x = 0). 
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